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Abstract 

The traditional thermodynamic Bethe ansatz (TBA) equations for the XXZ model at |A| > 1 
are derived within the quantum transfer matrix (QTM) method. This provides further evidence 
of the equivalence of both methods. Most importantly, we derive an integral equation for the 
free energy formulated for just one unknown function. This integral equation is different in 
physical and mathematical aspects from the established ones. The single integral equation is 
analytically continued to the regime |A| < 1. 

1 Introduction 

The thermodynamics of one- dimensional solvable models is generally determined by the solution 
to a set of socalled thermodynamic Bethe ansatz (TBA) equations Jl| . Some lattice spin models 
such as the XXZ chain, XYZ chain have been treated also by the quantum transfer matrix 
(QTM) method 0, [| |], ||, see also chapters 17 and 18 of [l|]. Correlated electron systems 
such as t-J model and Hubbard model, have also been treated by TBA and QTM methods 

The equations obtained by the QTM approach are quite different from those of TBA. 
However, the numerical results of the two methods for the free energies are the same. Mathe- 
matically the non-linear integral equations of [|lj and |J share similarities insofar as they can 
be interpreted as equations for dressed energies of elementary particles of magnon and spinon 
type, respectively. 

Recently, from the stand point of TBA one of the authors (MT) [l3j derived in the case of 
the XXZ chain a simple integral equation for just one unknown function. This integral equation 
is completely different in structure from those mentioned above. Here we aim at a derivation 
of this equation in the Quantum Transfer Matrix approach providing a more explicit as well as 
unified understanding of the structures and involved functions. 



To be definite, we first consider the region A > 1, 
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The thermodynamic Bethe ansatz equations for this model at temperature T are called Gaudin- 
Takahashi equations, [|TJ], [T5|] 
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The free energy per site is 
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From this equation MT |13[ derived 
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+ cot ^ - y + 2i] eM- 271 ^^ ^ ' 0])^^ (5) 
where the free energy is given by 

/ = -Tlnu(0). (6) 

The contour C is an arbitrary closed loop counterclockwise around where 2nQ,n ^ and 
±2z + 2nQ should lie outside of this loop. Furthermore this loop should not contain zeros of 
u(y). It is expected that u(y) has no zero in the region \Qy\ < 1. We show that these equations 
can be derived in the quantum transfer matrix approach. 



2 Quantum transfer matrix and fusion hierarchy models 

The quantum transfer matrix for this model is equivalent to that of the diagonal-to-diagonal 
transfer-matrix of the six-vertex model which is a staggered or inhomogeneous row-to-row 
transfer matrix, see below. The partition function Z = Trexp(— Tl/T) is given by 
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Then in the case N = 2M x integer, we have 

Z = TrT N , T((7i, cr 2 , <t 2 m; c 2 > — , o" 2 m) 

= A(a 1 a 2 ; a' 2M a' 1 )A(a 3 a i ; a' 2 a'^)...A[a 2M _ x a 2M \ <J 2M _ 2 (7 2M _ 1 ). 

The eigenvalue problem of this transfer matrix is a special case of the inhomogeneous six-vertex 
model on the square lattice. 

Consider an inhomogeneous six-vertex model with the following column dependent Boltz- 
mann weights: 

ai = pih(v + vi + 7)) 
k = piuj~ l \\{v + vi-i]) 
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Here L is the number of columns, h(w) is u, sin(-u) or sinh(w) depending on the anisotropy 
parameter. The transfer matrix T(v ) acts in a 2 L dimensional space, 
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The space is divided into subspaces characterised by the number of down spins k. Without 
loss of generality we can put k < L/2. In this subspace we can construct Bethe-ansatz wave 
functions with k parameters u±, Uk, 
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Imposing periodic boundary conditions the Bethe ansatz equations (BAE) take the form 
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The corresponding eigenvalue of the transfer matrix is given by 

Q(v) QO) 

Q(v) = Hh(v- Uj ). (12) 
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In order to solve the diagonal-to-diagonal transfer matrix we have to consider an inhomogeneous 
six-vertex model the Boltzmann weights of which are given by 



ai — Ci — 1, bi = b[ = for even I, 
ai = exp(— r f — ) sinh ( ) , bi = exp' 
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The conditions (P) are satisfied if we put 
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Putting r/ + V\ = 2ccm we have 

^ ) = ( h(2,)h(2 aM -2,) ) ' (16) 

The largest eigenvalue belongs to the k = M sector. The Bethe-ansatz equation for Uj, j = 
1, M are 
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The corresponding eigenvalue is given by 

T lW = ^(.-,)^ + ^ + ,)9^. (18) 

Due to the BAE fll7|) , the eigenvalue Ti(rr) is an entire function in the complex plane. The free 
energy per site is given by 

/ = — T lim lnTi(O). (19) 



The matrix T\(v) can be embedded into a more general family of matrices provided by the 
fusion hierarchy ||T7f , 



The eigenvalues Tj(v) as functions of v are all entire in the complex plane. It is easily seen 



that the following functional relations hold [17| 



Tj(v + ri)Tj(v -r}) = ip{v + (j + i)rj)<p(v - (j + 1)77) + T j+1 (v )Tj-i(v), 

T Q (v) = <p(v). (21) 

3 Derivation of Gaudin-Takahashi equation 

For A > 1 we put 
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h(u) = smu, = *0/2, = cosh" (—^Jj^fy), 

a>M = - tanh -1 f tanh 6 tanh - — ] , (22) 

In the limit of M — > oo we have 

= 0, Ma M = iJsinh0/(4T). (23) 

We transform the parameter v to x = iv/r). Then equations (|T6|) and ( p0|) turn into 
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T-(x) = y c -(^/T r /■ 2l)i) Q(x + (j + i)m(x-( J + i)i) 

These functions are all entire in the complex plane. Now we introduce a modified eigenvalue of 

%(x) = tM— ^M0)sinh0(l- MM ) sM 
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In contrast to the entire function Tj(x), Tj(x) has poles of order M at 

x = 2nQ + uui ± (1 + j — UAi)i- On the other hand, it has constant asymptotics 
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From (§1), we can find the following functional relation for Tj(x) 

Tj(x + i)Tj(x —i) — bj(x) + Tj_i(x)T i+1 (x), 

where we have defined 

-sin |(x + (j + 2uM)i) sin |(x — ji)\M 
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sin t(x + ji) sin t{x— (j - 2u M )i) 



Note that T Q (x) = 1 and bj(x), Tj(x) has poles at x = 2nQ + u M i ± (j — u M )i and 
x = 2nQ + umi ± ( j + 1 - um)% respectively. 
We define 



For these functions the following relations stand 
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As ^(x), j = 2, 3, ... has no pole or zero in —1 < Qx < 1, we find 

lnY^x) = s * (ln(l + Y^) + ln(l + Y j+1 )), j > 2. 
For Y"i(x) one must be careful that it has poles in —i, (1 — 2uu)i- Using 

T 2 (x + z)f 2 {x -i) = b 2 (x)(l + Y 2 {x)), 
and T2(x) has no zero or pole at — 1 < Qx < 1, we have 

lnf 2 (x) = s * (lnb 2 (x) + ln(l + Y 2 (x))). 
Using Yi(x) = T 2 (x)/bi(x) we have 

lnli(x) = — lnbi(x) + s * lnb 2 (x) + s * ln(l + Y 2 (x)). 
In the limit of M — > oo, the function bj(x) can be simplified 
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which has singularities at x = 2nQ ±ji. In the limit of M — > oo equations (|35"D , fl32|) , (|3lD are 
identical to (@). Substituting 



lnf = s * ln[(l + Yi(x)) /bi(x)] 



(37) 



into fll9|) we have Then the Gaudin-Takahashi equations are derived from the quantum 
transfer matrix method. (See also the treatment in [|T6| , |ITf for related models). 
Consider the M — > oo limit of the functions b j (xj . T j (xj as 

uAx) = lim T,(x), (38) 

Then from (|28"D we have the relation 

u\{x + i)u\(x — i) = bi(x) + v,2(x). (39) 

Note also the asymptotics Mi(±ioo) = 2 cosh h/T. We may assume the functions u\(x) and 
u 2 (x) have similar singularities at x = 2nQ ± 2i and x = 2nQ ± 3i, respectively. If we write 
© as 

■Ux(x + i) = bi(x)/ui(x — i) + m 2 (x)/'U 1 (x — i), (40) 

the LHS has singularities at x — i, —3i in the fundamental region (|5Rx| < Q). The first term 
of the RHS has singularities at x = i,—i,3i and the second term at x = 3i, — 3i — i. Then 



following the method in we get an integral equation for U\(x) 
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From the explicit expression of U\(x) (p6|), we see that the integral equation (|4]) is identical to 
the one obtained in flT3"|. The free energy is given by 
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4 Case A < 1 

In this case we have 
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In the limit of M — > oo we have 
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Putting x = iv/rj we obtain 
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Kuniba, Sakai and Suzuki jl7| succeeded in deriving the Takahashi-Suzuki equations for 
the thermodynamics of the XXZ model at h = 0, |A| < 1. The functions 
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are all periodic with periodicity 2poi. We have relations for Ti(x) and T 2 (x) 
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T\(x) satisfies 

Ti(±oo) = 2coshh/T. (49) 
By these two equations we can determine Ti(x) in the limit of M — > oo. In this limit &i(x) is 

h(x) = exp \- 2nJ ^ ne a 1 (x) ) , «i(x) = ^w^/^ (50) 
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We can assume that Ti(x) is expanded as follows 

T 1 (x) = 2coslA + YY- % - + fV 7 ^ -. (51) 
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Consider the contour integral around x = i giving the coefficients Cj 
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The first sum of the r.h.s. of (|51|) is 



J= i - ^ (a? - 2 ?W - 2%y Ti (y) 2vri J - x-y- 2np i - 2z Ti (y) 2m 

[9 9 27rJsin6> 1 dy 

= <t> — coth ~(x — y — 2i) exp — — cmy + z) . (53) 

7 2 2 V y 7 P[ Tfl vy ;J T 1 (y)27ri v 7 

The second sum is calculated in a similar way. Thus we find 
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+ coth -[x - y + 2i\ exp[ — — ai(y - i) ) ^^7^— , (54) 
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and the free energy is given by 

/ = -Tlnw(0). (55) 

Apparently these equations are analytical continuations of fl5|) and (§) if we replace by 26*. 
Then equation (|5|) treats the thermodynamics in a unified way 
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